On the discovery of doubly-magic 48 Ni 
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In an experiment at the SISSI/LISE3 facility of GANIL, we 
used the projectile fragmentation of a primary 58 Ni 26+ beam 
at 74.5 MeV/nucleon with an average current of 3 /iA on a 
natural nickel target to produce very neutron-deficient iso- 
topes. In a 10-day experiment, 28 7 42 Cr isotopes, 53 45 Fe iso- 
topes, 106 49 Ni isotopes, and 4 48 Ni isotopes were unambigu- 
ously identified. The doubly-magic nucleus 48 Ni, observed for 
the first time in the present experiment, is the most proton- 
rich isotope ever identified with an isospin projection T z — 
-4. At the same time, it is probably the last doubly-magic 
nucleus with "classical" shell closures accessible for present- 
day facilities. Its observation allows us to deduce a lower limit 
for the half-life of 48 Ni of 0.5 fis. 

PACS numbers: 21.10.Dr, 21.10.-n, 23.50.+Z, 25.70.Mn, 
27.40.+Z 

Doubly magic nuclei have a magic attraction for nu- 
clear physicists. These nuclei enabled them to interpret 
many properties of nuclei. They are spherical and ex- 
hibit a remarkable stability compared to neighboring iso- 
topes. However, one of the questions of current interest 
is whether or not the role of the magic numbers, well es- 
tablished along the valley of stability, remains important 
when an extreme excess of protons or neutrons is present 
in a nucleus. In fact, lighter naturally occuring isotopes 
have roughly the same number of protons and neutrons. 
Thus the lightest stable isotope of the element nickel is 
58 Ni with 28 protons and 30 neutrons, while 57 Ni is al- 
ready unstable and decays by (3 + emission. The excess 
of protons or neutrons in a nucleus is characterized by 
the isospin-projection quantum number T z which is half 
the difference between neutron and proton numbers. The 
most negative value for T z is expected to occur for 48 Ni. 

The recent observation of doubly-magic 100 Sn and 
78 Ni H showed that projectile fragmentation is most 
probably the route to 48 Ni and the question about its sta- 
bility attracted the attention of experimentalists as well 



as of theorists |4| |8j] . The special interest of this nucleus 
is that it is the only case of a doubly-magic nucleus in 
the entire chart of nuclei where the mirror nucleus, 48 Ca, 
is stable against particle emission, a fact which may al- 
low for interesting mirror symmetry studies. In addition, 
with 48 Ni, 56 Ni, and 78 Ni, nickel is most probably the 
only element possessing three doubly-magic nuclei. Fur- 
thermore, together with 45 Fe, 48 Ni is the prime candidate 
for the yet unobserved two-proton radioactivity, the emis- 
sion of 2 He from the ground state of an isotope If 
one excludes the super-heavy elements, where the magic 
numbers are not known experimentally, 48 Ni is also pre- 
dicted to be the last doubly-magic nucleus possessing 
shell closures known from spherical nuclei which is ac- 
cessible with present-day techniques. 

The search for doubly-magic 48 Ni started with the ad- 
vent of powerful projectile-fragmentation facilities. Thus, 
51 ' 52 Ni have been identified for the first time in 1986 Bfl 
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in an experiment at the LISE facility |10|. In this experi- 
ment, an average beam current of 50 nAcorresponding to 
1.5xl0 10 particles per second for the primary 58 Ni beam 
at 55 MeV/nucleon was available. 

The next two steps towards 48 Ni were taken at the 
fragment separator FRS |ll| at GSI. In a first experi- 
ment using the projectile fragmentation of a relativistic 
58 Ni beam at 650 MeV/nucleon (5xl0 7 particles per sec- 
ond) in 1993 [pi, 50 Ni was observed with 3 events. At 
the same time, extensive production cross-section mea- 
surements were performed which enabled us for the first 
time to reliably predict the production cross section of 
proton- rich isotopes from 58 Ni fragmentation (see below). 
After an increase of the primary-beam intensity of about 
a factor of 10 from SIS at GSI, a second experiment |l]J 
was carried out at the FRS using the projectile fragmen- 
tation of a primary 58 Ni beam at 600 MeV/nucleon to 
produce proton-rich isotopes. The experiment led to the 
discovery of 42 Cr, 45 Fe, and 49 Ni. 

Intensive developments on the ECR ion sources at 
GANIL made it possible to start a search for doubly- 
magic 48 Ni. In these developments, the metalic natural 
nickel used in the past was replaced by nickelocene, a 
chemical compound which allowed to treat nickel as a 
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gas in the ECR ion source. 

A primary beam of 58 Ni 26+ with an average current 
of 3 //A corresponding to 7.2 xlO 11 particles per second 
was accelerated by the GANIL cyclotrons to an energy 
of 74.5 MeV/nucleon. This beam impinged on a nat- 
ural nickel target of thickness 230.6 mg/cm 2 followed 
by a 2.7 mg/cm 2 carbon stripper foil, both installed 
on a rotating wheel between the two superconducting 
solenoids of the SISSI device. The secondary beams 
(Bpi = 1.7323 Tm) were transported through the Al- 
pha spectrometer to the LISE3 separator |Icfl . At the 
intermediate focal plane of LISE3, a shaped beryllium 
degrader (10.36 mg/cm 2 ) in conjuntion with the second 
dipole stage (B/?2 = 1-7000 Tm) allowed for a refined 
selection of the fragments of interest. 

At the first achromatic focal plane of LISE3, a micro- 
channel-plate (MCP) detector was mounted to yield a 
timing signal for each isotope. This detector had to with- 
stand a counting rate of about 5xl0 5 particles per sec- 
ond. It turned out that its efficiency was about 70%. The 
velocity filter at the end of the LISE3 device reduced the 
counting rate to roughly 500 particles per second. The se- 
lected isotopes were finally implanted in a silicon-detector 
telescope consisting of five detectors, the first three with 
a thickness of 300 each, the fourth with a thickness 
of 700 /jm and the last one with a thickness of 6 mm. 
The first detector (active area 600 mm 2 ) was used for an 
energy-loss and a time-of-flight measurement, whereas 
the second detector (quadratic shape with an area of 
3x3 cm 2 ) was position-sensitive (resistive read-out) and 
thus yielded an energy-loss and a two-dimensional po- 
sition measurement with a spatial resolution of about 
1 mm. The isotopes of interest were finally implanted in 
the third detector (600 mm 2 ) which measured the resid- 
ual energy and yielded another timing signal. Most of 
the lighter or less exotic particles were stopped in the 
first two silicon detectors. Thus the counting rate in the 
implantation detector was only about 10-15 particles per 
second. The last two detectors served as veto detectors 
and had a size of 15x40 mm 2 and of 600 mm 2 . 

This detection set-up enabled us to measure 10 inde- 
pendent parameters to identify the isotopes for each im- 
plantation event: three energy signals in the first three 
silicon detectors, two veto signals in the fourth and the 
fifth detector, positions in x and y, and three independent 
flight times (TOF). The TOFs were measured between 
the MCP detector and the first silicon detector, the first 
silicon detector and the cyclotron radiofrequency, as well 
as between the third silicon detector and the radiofre- 
quency. 

To analyse the data, the average value for each of these 
10 parameters as well as its full width at half maximum 
(FWHM) was determined for each isotope based on the 
experimental observation. In order to calculate the pa- 
rameters for particle-unstable isotopes not present in the 
spectra, interpolations and extrapolations from the mea- 
sured values of neighboring nuclei were used. The ex- 
pected average value for isotopes with isospin projections 



from T z — -2 to T z — -9/2 and with nuclear charges from 
Z = 22 to Z = 29 was calculated for each parameter. 
To accept an event, all parameters had to lie within a 
window of two FWHM. However, for the graphical rep- 
resentation, these cuts were only applied for parameters 
not represented in the spectrum. Figure 1 shows the re- 
sult of this analysis as a plot of the energy loss in the 
first silicon detector as a function of the TOF between 
the MCP detector and the first silicon detector. Due to 
the efficiency of the MCP detector, part of the statistics 
is lost in this spectrum. Nevertheless, two events of 48 Ni 
are clearly observed. In particular, no background events 
are observed in the whole region analysed between the 
unbound Cu at the upper left end and the unbound 
35 Ti at the lower right end. In addition, 77 events of 
49 Ni, 29 events of 45 Fe, and 164 events of 42 Cr fullfil 
the stringent conditions imposed on this spectrum. Al- 
ready from this spectrum we conclude therefore that the 
new doubly-magic nucleus 48 Ni has been unambiguously 
identified. 
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FIG. 1. Two-dimensional identification plot of energy loss 
in the first silicon detector versus the TOF between the MCP 
detector and the silicon detector. To produce this spectrum, 
software conditions were applied to all parameters except 
those represented here (see text for details). Due to the effi- 
ciency of the MCP detector, part of the events are lost in this 
spectrum. 

To circumvent the lack of efficiency of the MCP detec- 
tor, we analysed the data also as a function of the energy 
loss in the first silicon detector and the total energy as 
the sum of the signals in the first three silicon detectors. 
For this analysis, we used the information from the MCP 
detector only, if it was present, i.e. if the TOF signal 
of this detector and the first silicon detector was above 
zero. The result of this analysis is shown in figure 2. In 
this spectrum, we observe four counts which we can un- 
ambiguously attribute to Ni. The number of counts for 
49 Ni, 45 Fe, and 42 Cr are 89, 48, and 264, respectively. 

This spectrum, however, demonstrates also the impor- 
tance of the MCP detector, as all background counts are 
events without MCP timing. In fact, we identify one 
background count at the positions of the unbound 51 Cu, 
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45 Mn, 42 V, and 41 V, whereas the other counts could be 
rejected due to their wrong energy-loss/energy relation 
in this spectrum. However, none of these background 
counts are accumulated in a particular region in this 
energy-loss/energy matrix, like it is the case for the four 
48 Ni events. Therefore, it seems rather unlikely that more 
than one of these four events is background. 
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FIG. 2. Two-dimensional identification plot of energy loss 
in the first silicon detector versus the total energy calculated 
as the sum energy of the first three silicon detectors. As for 
figure 1, software conditions were applied to all parameters 
except those represented here (see text for details). Less ex- 
otic isotopes were cut either by software conditions or by the 
transmission of LISE3. 

From figures 1 and 2, we conclude that the new doubly- 
magic nucleus 48 Ni has been observed for the first time. 
The observation of 49 Ni, 45 Fe, and 42 Cr confirms the re- 
sults from the 1996 GSI experiment [|l3| where these iso- 
topes were observed for the first time with much lower 
counting rates. 

The number of counts of the present experiment for 
the four most proton-rich isotopes 48 Ni, 49 Ni, 45 Fe, and 
42 Cr is given in table 1. In order to have all counts for 
each nucleus to calculate the production cross sections, 
we used three-FWHM windows to determine the number 
of counts. To estimate their production cross section we 
calculated the Alpha+LISE3 transmission for these iso- 
topes with three different simulation codes, the programs 
LISE @, LIESCHEN @, and INTENSITY @. The 
INTENSITY code has no feature to treat the velocity fil- 
ter of the LISE3 separator. However, it turned out in the 
simulations with the two other codes that the velocity- 
filter transmission for the four isotopes discussed here is 
100% and thus INTENSITY calculations yield the total 
transmission in these special cases. In table 1, we give 
the average transmission from the three simulations. The 
error is estimated from the difference of the predictions. 
Except for 42 Cr which is close to the acceptance limit of 
LISE3 at the intermediate focal plane, the different codes 
predict transmissions within a window of about 50%. 

The production cross sections for 48 Ni, 49 Ni, 45 Fe, and 



42 Cr in the present experiment can be deduced from this 
information, the ™°*Ni target thickness of 230.6 mg/cm 2 
and a total beam dose of 4.2xl0 17 particles. These cross 
sections are presented in table 1. They are compared to 
the new EPAX formula pj], an empirical parametriza- 
tion of projectile fragmentation cross sections. We find 
nice agreement between experimentally determined cross 
sections and predictions from EPAX, especially for 48 Ni 
and Ni. However, the experimental production cross 
section for 45 Fe is more than a factor of two lower than ex- 
pected. From the cross-section systematics measured for 
the fragmentation of a 58 Ni primary beam |12[ , a decrease 
of a factor of two is not expected between ^Ni and 45 Fe. 
A possible explanation could be that part of the 45 Fe ac- 
tivity decays in-flight with a half-life of roughly 0.5 /is. 
However, we can not exclude that the production cross 
section for 45 Fe might be smaller than expected or that 
the transmission calculations have a larger error. There- 
fore, additional information, namely measurements of its 
decay as well as a measurement similar to the present 
experiment, but with a significantly shorter flight time, 
are necessary to confirm or reject this hypothesis. 

48 Ni is predicted particle unstable by all commonly 
used mass models. As it is an even-Z nucleus, these mod- 
els foresee that the one-proton emission is energetically 
forbidden and thus 48 Ni may decay by two-proton emis- 
sion. However, the model predictions for the Q value 
for two- proton emission, Q2p, vary between less than 
1.0 MeV and more than 3.0 MeV (see figure 3). Us- 
ing barrier-penetration calculations this gap yields par- 
tial half- lives for two-proton decay between 10~ 14 s and 
Is. 





48 Ni 


49 Ni 


45 Fe 


42 Cr 


number of counts 
transmission (%) 
cross section (pb) 
EPAX (pb) 


4±2 
9.8±1.2 
0.05±0.02 
0.06 


106±22 
6.9±1.6 
1.8±0.6 
1.5 


53±9 
6.7±0.3 
0.9±0.2 
1.8 


287±30 
0.7±0.6 

47±i 5 
69 



TABLE I. Pertinent information concerning the produc- 
tion cross sections of 48 Ni, 49 Ni, 45 Fe, and 42 Cr. The first 
line gives the number of counts experimentally observed. The 
numbers given correspond to 3-FWHM windows for each pa- 
rameter (see text). The error bars are the sum of statistical 
and systematic errors. The systematic errors have been deter- 
mined by varying the conditions used to identify the isotopes. 
The transmissions in the second line are averages from three 
different model calculations with LISE @, LIESCHEN jfi), 
and INTENSITY (le) . The experimental cross sections (third 
line) are determined with a primary beam dose of 4.2xl0 17 
particles and a nat Ni target thickness of 230 mg/cm 2 . In ad- 
dition, a dead time correction of 13% is applied. Losses due 
to secondary reactions in target and degrader are below the 
one percent level and are not taken into account. The last 
line give s the predictions from the new EPAX parametriza- 
tion 
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However, a closer inspection of the mass models used 
in figure 3 shows that "local models" , i.e. models which 
either have been tailored to calculate masses in the region 
around 48 Ni ||-[7) or which use local mass relations like 
the isobaric multiplet mass equation Jl8fl or the Garvey- 
Kelson mass relation jl9| , predict Qi v values which clus- 
ter at about 1.3 MeV. Some of the "global models" which 
consistently predict masses for a large part of the table 
of isotopes give also masses in this Q-value range [ pJ[ . 
However, most of them lie above 1.5 MeV and predict 
therefore much shorther half-lives. In microscopic calcu- 
lations like those of Nazarewicz and co-workers M , the 
masses depend on the Skyrme force used. 
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FIG. 3. Barrier-penetration half-life as a function of the 
two-proton Q value, Q2 P , for 48 Ni. The barrier penetration 
was calculated by assuming a spectroscopic factor of unity. 
Different model predictions |p|-PJl^-p^ were used for Q2 P - 
The experimental observation of ^i implies a half-life longer 
than 0.5 /ts yielding an upper limit for Q2 P of about 1.5 MeV. 

From the experimental observation of Ni, we can de- 
duce a limit of its half-life. Assuming that 49 Ni does 
not decay in flight, we can use the number of counts 
observed for this isotope to estimate the number of 
counts expected for 48 Ni. The transmissions for 48 > 49 Ni 
(see table 1) and a cross-section decrease of a factor of 
26.4±4.7 pj] allow us to deduce a number of counts ex- 
pected for^ 8 Ni of 5.7±2.2. With these figures and a 
total flight time of 1.32 /is, we obtain a half- life for 48 Ni 
of (2.6j]^ c 1 )/is. We determine thus a lower limit at the la 
level of 0.5 /is for the life time of 48 Ni. This value trans- 
lates into an upper limit for Q2 P of about 1.5 MeV if we 
assume a spectroscopic factor of unity. Combined with 
the mass excess of (0.755±0.354) MeV for 46 Fe |23), we 
determine therefore an upper limit for the mass excess 
of 48 Ni of 17.2 MeV. One should keep in mind that this 
determination is based on an extrapolated mass excess 
for 46 Fe and a spectroscopic factor of unity which make 
it only a rough estimate. 

The /3-decay half-life of 48 Ni is predicted to lie within a 
range of 5.9 ms and 11.3 ms ^j^4[|. Therefore, if the Q2 P 
value is, as predicted by several mass models, of the order 
of 1.3 MeV, which yields a barrier-penetration half- life of 



about 10 _5 s, 48 Ni should decay by two-proton ground- 
state emission. Such a half-life should be accessible in 
future experiments. 

In summary, we observed for the first time 4 events of 
the doubly-magic nucleus Ni. The experimental obser- 
vation implies that its half-life is longer than about 0.5 /is. 
Using a spectroscopic factor of unity, this translates into 
an upper limit for Qi v of 1.5 MeV. Several global mass 
models are in disagreement with this value, whereas lo- 
cal mass relations are in nice agreement with the present 
observations. Future measurements of the decay of 48 Ni 
which should be feasible after additional improvements 
of the primary-beam intensity should yield more detailed 
information on its mass and its decay properties. 
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